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The linear response of a superfluid, rotating uniformly in a cylindrical container and 
threaded with a large number of vortex lines, to an impulsive increase in the angular 
velocity of the container is investigated. At zero temperature and with perfect pinning of 
vortices to the top and bottom of the container, we demonstrate that the system oscillates 
persistently with a frequency proportional to the vortex line tension parameter to the 
quarter power. This low-frequency mode is generated by a secondary flow analogous 
to classical Ekman pumping that is periodically reversed by the vortex tension in the 
boundary layers. We compare analytic solutions to the two-fluid equations of |Chandler| 
& Baym (1986) with the spin-up experiments of Tsakadze & Tsakadze ( 1980| ) in helium 
II and find the frequency agrees within a factor of four, although the experiment is not 
perfectly suited to the application of the linear theory. We argue that this oscillatory 
Ekman pumping mode, and not Tkachenko modes provide a natural explanation for 
the observed oscillation. In neutron stars, the oscillation period depends on the pinning 
interaction between neutron vortices and flux tubes in the outer core. Using a simplified 
pinning model, we demonstrate that strong pinning can accommodate modes with periods 
of days to years, which are only weakly damped by mutual friction over longer timescales. 
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1. Introduction 

The linear response of a rapidly rotating Navier-Stokes fluid to an impulsive increase in 
angular velocity of its container, or ‘spin-up’ has been extensively studied. In the seminal 


work of Greenspan & Howard (1963), three phases were identified: formation of a viscous 


boundary layer, Ekman pumping, and the damping of residual inertial oscillations by 
viscosity. Co-rotation between the interior fluid and container is established by Ekman 
pumping, during which a secondary flow recycles fluid from the boundary layers into the 
interior. The time-scale for the spin-up, known as the Ekman time, is proportional to the 
viscosity to the minus half power, and achieves co-rotation much faster than viscous dif¬ 


fusion. Since the work of Greenspan & Howard (1963), spin-up has been studied in fluids 
with density stratification ( Pedlosky|1967 Walin|1969 l, stratified and compressible fluids 
(Abney & Epstein 1996; van Eysden fe Melatos||2008 ) magnetized plasmas (Loperfl97T 
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Easson||1979 van Eysden|2014), and multi-component fluids (|Ungarish||1990 Amberg & 


Ungarish|l993) and different geometries (Clark et al||1971 van Eysden fe Melatos||20T3 


2014). For a review of spin-up in classical fluids, the reader is referred to (Benton & Clark 
|I974 ). 

The spin-up of superfluids like helium II has also received considerable attention both 


experimentally (Tsakadze & Tsakadze 1980) and theoretically (Campbell & Krasnov 
1982 Adams et al.||198l) Peradzynski et aZ.||1990 ) with the aim of shedding light on the 


interaction between the quantized vortices and the container walls. A rigorous treatment 


was provided by Reisenegger (1993), who solved the two-fluid equations of Chandler & 


Baym (1986) between slowly accelerating parallel plates, and showed that the Ekman 


time is reduced by a factor depending on the superfluid mutual friction coefficients and 
the normal fluid density fraction. This was generalized to the impulsive acceleration 


of containers of arbitrary shape ( 

van Eysden & Melatos 

2013) 

and to include the self- 

consistent response of the container I 

van Eysden & Melatos 2014 

to facilitate comparison 

with helium II experiments ( 

van Eysden & Melatos 

2011 

201 

2). However, theoretical 


work to date on the spin-up of two-component superfluids such as helium II assumes 
smooth-walled containers and neglects the effects of pinning. 

The spin up of helium II at zero temperature, where the normal fluid component 
vanishes and cannot facilitate Ekman pumping was also investigated by |Reisenegger| 
(1993). By introducing a frictional force to account for the sliding of vortex lines at the 
boundary, Reisenegger (1993) showed that the superfluid spins up via an Ekman-like 
secondary flow with a timescale that depends on the strength of the frictional force. 
However, only the response of a single-component fluid to slowly accelerating parallel 
plates was studied. 

There is good reason to expect that the impulsive acceleration problem produces dif¬ 
ferent physics than that of the slowly accelerating parallel plates. When an ideal magne- 
tohydrodynamic plasma is slowly accelerated (i.e., the time-scale for acceleration of the 
container is the slower than the Alfven crossing time), the plasma moves with the con¬ 
tainer as a rigid body. However, in response to an impulsive acceleration, Alfven waves 


are excited that produce persistent oscillations of the container and the plasma (van Eys- 
denj.2014). A final state of co-rotation between the container and plasma is inconsistent 


with energy conservation if the fluid is dissipation-less, hence the system oscillates persis¬ 
tently. Similar arguments apply for a superfluid at zero temperature rotating uniformly 
with a high density of vortex lines. 

Two primary motivations exist for this study. The first is a series of experiments 
performed by Tsakadze & Tsakadze (1980), in which the angular velocity of containers 
filled with uniformly rotating helium H was impulsively increased. The subsequent motion 
of the container, responding freely to the hydrodynamic torque of the contained fluid, was 
then recorded. In experiments where the container was coated with powder to facilitate 
vortex pinning, a sinusoidal oscillation of the container was observed. These oscillations 
are have been interpreted as Tkachenko modes, however the dependence of the oscillation 
period on vessel radius predicted for these modes is inconsistent with measurements, 
whereas the columnar nature of the Ekman pumping mechanism presented here is not. 

The second application is the recovery of pulsar glitches; the original motivation for 
experiments of Tsakadze & Tsakadze (1980). Glitches are tiny, impulsive increases in the 
rotation frequency of the neutron star crust, typically followed by a quasi-exponential 
recovery that is believed to be associated with the response of the interior superfluid 


Baym et al. (1969a). Although the superfluid neutron vortices are believed to pin to 


lattice sites in the crust and to flux tubes arising from type II superconductivity of 
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protons in the outer core, the Ekman pumping mechanism identified here may give rise 
to long-period oscillations if vortex creep is active. 

In this paper, we revisit the problem of superfluid spin-up, by solving for the response 
of a uniformly rotating superfluid to an impulsive acceleration of the container. We focus 
on cylindrical geometry, where the vortices are strongly pinned to the top and bottom 
of the container. We apply the traditional Laplace transform techniques of Greenspan 


& Howard (1963), but solve self-consistently for the motion of the container as in 


van 


Eysden & Melatos (20141 and van Eysden (2014). We demonstrate the presence of a 


low-frequency oscillation mode, with a period that scales as the vortex line parameter to 
the one-fourth power. This mode arises in the fast rotation limit, analogous to classical 
Ekman pumping. A secondary flow is present, which is periodically reversed by the 
tension in the vortex lines as they are sheared in the boundary layer. By solving the two- 
fluid equations of Chandler & Baym (1986), we predict an oscillation period of 10 seconds 
in helium II, compared with that of 40 seconds observed in the experiments of |Tsakadze| 


& Tsakadze (1980). A direct comparison cannot be made because the experiment is not 


perfectly suited for the application of the linear theory of Chandler & Baym (1986). Using 
a simplified pinning model, we predict that oscillations in neutron stars can have periods 
of days to years if the pinning of neutron vortices to flux tubes in the core is strong. 
These oscillations are weakly damped and expected to be observable over timescales 
much longer than the period. 

The paper is structured as follows. In f|2] we linearize the equations of Chandler fe 
Baym (1986) using a common Ekman pumping ansatz. These equations are solved in f|3 
for a superfluid at T = 0 and perfect vortex pinning of vortices at the top and bottom 
of a cylindrical container, where it is shown that the solution is oscillatory and has 
an Eknran-like secondary circulation. In ; |T[ the full two-fluid equations are solved and 
compared with the experiments of Tsakadze fe Tsakadze| (1980) in helium II. Neutron star 
applications are considered in (j5j In (J 6 ]we summarize our conclusions. Technical details 
of the solution are presented in the appendices along with some important limiting cases. 


2. Governing equations 

A convenient description of superfluids such as helium II, rotating with a high density 
of vortex lines is given by the Hall-Vinen-Bekharevich-Khalatnikov (HVBK) equations. 
The fluid comprises two components: a ‘normal’ component, denoted by subscript n, 
with viscosity 77 , and an inviscid ‘superfluid’ component, denoted by subscript s.[f]Under 
rotation the inviscid component forms a dense array of quantized vortices, which are 
smooth-averaged in the hydrodynamic approximation, endowing the inviscid component 
with a macroscopic vorticity. The vortices mediate interactions between the normal and 
superfluid components, giving rise to a mutual friction force. An informative introduction 


to the HVBK equations is provided in Henderson & Barenghi (2000). A more general set 
of equations, derived by Baym & Chandler ( |1983 ); Chandler & Baym (1986), includes 
restoring forces experienced by the vortex lattice when vortices are displaced from their 
equilibrium configuration. This force produces Tkachenko oscillations, which are expected 
to be observed in superfluid experiments. Assuming that the fluid is incompressible, the 


equations of Chandler & Baym (1986), written in HVBK form and in the laboratory 


I In a neutron star, the ‘normal’ and ‘superfluid’ components refer to a proton-electron plasma 
and neutron superfluid, respectively. This is discussed further in Sj5] 
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frame are 


dv n 

dt 

dv s 

dt 


V n ■ VVn 

- v s ■ Vv s 

V • v n 

V • v s 


— + —F + v n S7 2 v n , 

p 

= -Vps - i - a - —F , 

P 

= 0, 

= 0, 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 

(2.4) 


where u rajS , p niS and p n ^ s are the macroscopic velocities, densities and reduced pressures 
of the normal and superfluid components, respectively. Throughout our analysis, both 
the normal and superfluid components are considered incompressible; the validity of this 
assumption is assessed in f|4] and f|5] for helium II and neutron stars respectively. The 
mutual friction force is 

F = x [w, x {v„ — v s ) — t — cr] 

+ ^B'[u) s x(v n -v s )-t-(r], (2.5) 

and 

u) s = V x v s , (2.6) 

is the macroscopic vorticity of the superfluid and uj s = uj s /\ijO s \ is the vortex line direction 
vector, and B and B' are dimensionless mutual friction coefficients. The vortex tension 
force per unit mass is 

t = v s uj s x (V x u) s ) . (2-7) 

where the vortex line tension parameter is given by 

■'• = £ iog O' (2 - 8 > 


Here T = hn/m is the quantum of circulation, m is the mass of one helium atom, b 0 is 
the inter-vortex spacing and ao is the size of the vortex core. The kinematic viscosity 
v n is defined as the shear viscosity divided by the normal fluid, density, rj/p n . The final 
parameter er in (2.1)—(2.41 comes from the theory of Baym & Chandler (1983) and 


describes the restoring force of the vortex lattice in response to shear deformations. It 
has the form 

h\u3 s 


8m 


[2V ± (V -0 - Vi£] , 


(2.9) 


where £ is the vortex line displacement vector and Vj_ is the two-dimensional gradient 
operator. Both £ and Vj_ are two-dimensional in the plane orthogonal to the angular 
velocity of the background superfluid flow, i.e., uj s ■ £ = ui s ■ = 0. Equation (2.9) 


only applies to linear deformations of a rectilinear vortex array, hence (2.1)—(2.4) are 
only valid in the linear approximation when er is included. When er = 0, (2.1)—(2.4) are 


the HVBK equations, which describe non-linear flows including quasi-classical turbulence 
( Henderson et «L~|| 1995 Peralta et a(.||2008 ). 

The vortex lines obey the vorticity conservation law 


dU)s Y7 ( \ 

—— = V x (v L x (jj a ) 
ot 


( 2 . 10 ) 


where in the linear approximation the perturbation to the vortex line velocity is given by 


d^/dt. Taking the curl of (2.21 and comparing the result with (2.101 gives the equation 
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of motion for the vortex lines 

u) s x (v L - v s ) = t + cr + —F . 

P 


Equations (2.1) and (2.2) can be combined into an equation for the total fluid, 
d 

{Pn^n H - Ps'^s') ~\j {Pn'^ni'^nj H“ Ps'^si'^sj^) = ^ j^ij i 

where 

i\, - -p<\, + rr ■ i . 

The contributions to the stress are the total pressure 

P = PnPn + PsPs - p s V s \u>s\ . 

the viscous stress 


Tij = PnVn (Vi 


iy + Vj 'i 


i), 


the vortex line tension, 


Ty j — Ps^S |^-s {0 s iCj S j Sij) , 


( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 

(2.14) 

(2.15) 

(2.16) 


and the stress arising from the displacement of vortices from the equilibrium configuration 
in the lattice, 


rjit _ 

1 ij — 


Psh\u3 s 

8 TO 


[V_u£j + Vj — 36 ij (dkfk)] 


(2.17) 


The term (2.17) is responsible for Tkachenko oscillations. 


To study the coupled response of a superfluid and its container, we consider two infinite 
parallel plates with separation 2 L. This geometry has a long history in the study of the 


spin up of rapidly rotating fluids in geophysics (Greenspan & Howard 1963 Pedlosky 


1967 


Walin|1969), magnetized plasmas (Loper|1971 Easson fe Pethick|1979 van Eysden 


2014) , condensed matter ( Reisenegger||1993 ) and astrophysics ( Abney fe Epstein||1996 ). 
At times t < 0, the superfluid and its container rotate rigidly and uniformly about the 
cylindrical axis with angular velocity Q. At time t = 0, the magnitude of the angular 
velocity of the container is impulsively increased to fl(l + e), where the e 1 is the 
Rossby number. For t > 0 the container and fluid are left to evolve freely, and we solve 
self-consistently for the coupled motion of the container and plasma as in previous studies 


(van Eysden & Melatos 

2013 

to 

o 

I — 1 

van Eysden 

2014 


and neglected, however they typically play a secondary role in spin up problems like that 
considered here. Side wall effects in helium II experiments are discussed in §4.3| 

To solve for the motion of the container, we invoke angular momentum conservation 
between the container and fluid 

Ic-jf = - <b XX (h- T) dS + T ext , (2.18) 

where I c is the moment of inertia of the container, n is the unit vector normal to the 
boundary and dS is an element of area on the boundary. The first term on the right hand 


side of (2.18) is the hydrodynamic torque exerted on the crust by the fluid arising from 


viscous stresses and stress exerted by the vortex array, where T is given by (2.12). The 


second term is an external torque which may arise from, e.g., friction in the apparatus 
for superfluid experiments or the magnetic dipole torque in pulsars. 

The normal fluid co-rotates with the container, giving the boundary condition 


»„=9,xa). 


(2.19) 
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where x is the radial vector and f 2 c (f) is the angular velocity of the container, which is 
a function of time. Equation (2.19) embodies the usual no-slip boundary conditions for 
viscous flows. Following Reisenegger (1993), for the superfluid we choose the following 
boundary conditions 


n x [p s \u: s \L'y (v L - fl c x x) ± (n ■ T)] = 0, 


fi ■ v s = 0 , 


( 2 . 20 ) 


on surfaces that intersect vortex lines, where fi is the unit vector normal to the surface. 
The dimensionless constant 7 governs the rate of vortex creep at the boundary. In the 
limit 7 —► 00 , the vortices are pinned to the boundary and when 7 = 0 , the vortices exert 
no stress on the boundary, i.e., they are freely sliding. For the normal component of the 
superfluid we require no-penetration. 

In cylindrical coordinates (r,<f>,z), the initial conditions are 


v n ,s( 0 ) = rSl$, 


J 1 C ( 0 ) = fi(l + e)z , 


( 2 . 21 ) 


i.e., we assume the two fluids are initially co-rotating. Strictly speaking, the initial ve¬ 


locity for the fluid components obey ( 2 . 21 ) everywhere except in an infinitely thin region 


adjacent to the boundary where it is spun up by the container. 


For e <C 1, equations (2.1)—(2.21) can be linearized by perturbing around an equilib¬ 


rium rotating with uniform angular velocity fl about the z-axis. The external torque is 
also taken to be aligned with the rotation axis, i.e., T ext = T ex tZ. The geometry and 
initial conditions are axisymmetric, and the resulting flow axisymmetric. The following 
substitutions are made for the velocity and pressure fields 


v ntS (r, z , t) —> rVLcj) + eflL 
vl(t, z , t) —► r£l(j> + ef ILr* 


r* d * f + r*V n , s (j) - 2xn, s z 


dU £ , dV( 

—-r -I- 

dt dt 


p n Ar,z,t)^^+en 2 L 2 ( r -^^+2Q n , s ) , 


( 2 . 22 ) 


where Xn,s( z *,t*), V n , tS (z*, t*), Q n ,s(z*,t*) and P n , s (t*) are all dimensionless quantities. 
The functions r\n,s are stream-functions for the secondary flow. The asterisked quantities 
are defined as r* = r/L , z* = z/L , t* = ilt. Equation (2.22) is essentially the “von- 


Karman similarity” form and is typically used in employed in studies of spin-up between 
parallel plates ( Greenspan fc Howard||1963 Easson|[l979 Reisenegger||1993 van Eysden 
Sz Melatos|2013 ). The ansatz (2.22| automatically satisfies the continuity equations (2.3) 
and (2.4) and the conditions for rotational equilibrium for the background flow. The radial 


dependence of the azimuthal velocity is motivated by the boundary conditions (2.19), 


and chosen for the other variables such that r vanishes from the resulting equations in 
the most general way. Under these assumptions, the only non-vanishing component of 
the external torque is in the z direction, hence 


n c (t*) —► Viz + ef, 


(2.23) 


where / is a function of t* only. Henceforth, we drop the asterisk notation so that all 
variables are dimensionless. 


Substituting (2.22) into the normal fluid momentum equation (2.1), we obtain 

d 2 pJB\ d'Xr, 

P 


° [dt E dz 2 


_ ( 2 - PsE ) V 
dz \ P ' 
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p s B 

^P~ 

d 


2 -E, 


d 2 


° [dt E dz 2 


p s B 

' 2 p 


dz 2 

d 2 , 

P 
d 2 


2-E s 




<9z 2 

Xn 


dxs 

cb 

K 

K + 

<9Qn 

dz 


p s B' 
2 P 

2 - 


2-E, 


d 2 


dz 2 

p s -B'\ 


Vi + P„ 


p s 5' 

2 p 


P 

2 — E, 


dz 

d 2 \ d Xs 


dz 2 / 9^ 


7 

(2.24) 

(2.25) 

(2.26) 


and from (2.2) we have for the superfluid 

d 2 


o = 


o = 


PnB 


dt 2 p 

PnB dXn 

p dz 
d_ p n B 
dt 2 p 

PnB 


2-Es 


dz 2 


- 1 - 


PnB' 


V n 


2-Es 


cf_ 

dz 2 


dXs 

dz 


Vs + 1 - 


PnB' 

2 P 


2-Es 


dz 2 


Vs 


PnB’ 

2 P 


2-Es 


d 2 \ d Xs 


dz 2 J dz 


-V n 


p n B' dXn 


P 


0 = 


dXs 


p dz ’ 

dQ s 

dt dz ’ 
where we define the dimensionless parameters 

„ Vn 

E = 


L 2 Q 


Es = 


L 2 n 


(2.27) 

(2.28) 

(2.29) 

(2.30) 


The dimensionless parameter E is the Ekman number and is a ratio of viscous forces and 
the rotational inertia in the flow. By analogy, we define the superfluid Ekman number 
E s , which is a ratio of the vortex line tension and rotational inertia in the flow. Equation 
(2.18) becomes 


df K 

^7=4 — 
dt p 


PsE. 


d 2 Xs 

dz 2 


PnE 


VJVn 

dz 


at z = ±1, where we define 


I< = 


(1 + K) , 

npR 4 L 


a = 


'Text 

(El 2 Itc 


(2.31) 


(2.32) 


In (2.31), the dimensionless parameter K denotes the ratio of the moments of inertia of 
the fluid and container, and a is the dimensionless external torque where I to t = ^c(l + i4) 
is the total moment of inertia of the fluid and container. 


The boundary conditions for the normal fluid (2.19) become 


V n -rf = 0, 
dXn = 
dz 

Xn = 0, 


(2.33) 

(2.34) 

(2.35) 


at z = ±1. For the superfluid, using the linearized forms of (2.11) and (2.2) to eliminate 
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the vortex line velocity and mutual friction, ( 2.201 becomes 


d dxs 
dt dz 


dV 1 ± E 1 dV 1 

dt 7 dz 


2 f± 


E s d 2 \s 
7 dz 2 


+ Pa 


Xs 


at z = ±1. The initial conditions (2.21) become 


V n (z, 0) =0, 
V a (z,0) = 0, 

m = i • 


o, 

o, 

o, 


(2.36) 

(2.37) 

(2.38) 

(2.39) 

(2.40) 

(2.41) 


The vortex line displacements, £, can be calculated by substituting (2.22) into the vortex 
line equation of motion (2.11) and using (2.2). We find 


2 d^ = d_dx^ 
dt dt dz 
dUs _ dV s 
dt dt ' 


(2.42) 

(2.43) 


For the initial conditions we require 


Ut(z, 0)=0, (2.44) 

V((z, 0) = 0. (2.45) 


3. Pure superfluid 

3.1. Exact solution 

To illustrate the Ekrnan pumping mechanism which is the principal result of this paper, 
we consider a superfluid at T = 0, where p n = 0. To solve the initial value problem, we 
take the Laplace transform, 


X(z, s) = / X(z, t)e 


*dt. 


(3.1) 


Equations (2.27) and (2.28) are 

d\s 
S dz ' 


d 2 \ ~ 

E sg^) V s + P s = 0, 

(3.2) 


(3.3) 


Equation (2.29) is not required to solve the system, and is only necessary if one desires 

(3.4) 


to calculate Q s for completeness. Taking the derivative of (3.2) we obtain 

d 2 Xs 


dz 2 


- 2 -E, 


d 2 \ d \4 


dz 2 1 dz 


= 0. 


Equations (3.3) and (3.4) can be solved for \s and V s ; P s is then obtained from (3.2). 
The equation for the container, (2.31) becomes 


sf - 1 = T KE s 


d 2 Xs 

dz 2 


(3.5) 
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at z = ± 1 , and the boundary conditions are 


7 dz 


sK±—^ = 0, 


s fi-2/±&f%+P» = 0 , 

oz 7 

= 0 , 


at ^ = ±1. The solution to (3.3)—(3.8) is 

_ 2 if 

V s = {C .|_ [«_s(coshK_ — coshK_^) + (E s /^)k 2 _ sinh«;_] , 

— C_ [k+s (cos1ik + — cosh/c + z) + (E s / j)k+ sinhfc + ] } , 

2 f 

Xs = -r— [C+ (sinh — z sinh «;_) + C_ (sinh n + z — z sinh re + )l , 
As 

P s = -^4- (4 + s 2 ) ((7+ sinh «_ + C_ sinh k + ) , 

As^ 

f = A 

J As ’ 

where 


k± = 


2 ± is 


E„ 


C± = k± cosh k± + 


(E s k± =f 217 ) sinh«± 


7 s 


A = C+ 
A = C+ 
+ C _ 


C_ 


c_ 


+ c_ 


2 E s k 2 _ (1 + I \) sinh «_ 


C+ + 


2E s k 2 + sinh re + 


2E s k 2 + (1 + K) sinh n + 

C + H-7 

s z 


(3.6) 

(3.7) 

(3.8) 


(3.9) 

(3.10) 

(3.11) 

(3.12) 

(3.13) 

(3.14) 

(3.15) 

(3.16) 


To obtain the response of the container, we must find the inverse Laplace transform of 
/. This is readily done by realizing that there are simple poles at s = 0 and at the zeroes 
of A. The result is 


/(*) = 


where s n are the roots of A and 


J s =0 


R(s) = 


7 ~^R(s n )e s 


A(s) 

„dA 
* ds 


(3.17) 


(3.18) 


The sum in (3.17) implies summing over all the zeroes of A.[f]ln general, the eigenvalues 
occur in conjugate pairs. For 7 -> 00 we have Re(s n ) = 0 and Im[i?(s„)] = 0, and 
the modes in (3.17) are purely oscillatory. When 7 is finite the friction force results in 
dissipation and s n and R(s n ) have real and imaginary components. 


f The solution to the initial value problem and could also be obtained from a mode expansion 
and appropriate orthogonality relations, however, the Laplace transform more readily gives the 
result. 
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Figure 1. Dimensionless angular velocity of the container, f(t), for E s =0.1 (blue curve) and 
E a = 0.01 (purple curve). We take K = 1 and 7 —► 00 . 



100 




0.25 



90 


• 





80 




0.2 

• 


70 


• 





60 


• 


0.15 


s 

50 






40 


• 

u oj 

0.1 



30 


• 

& 




20 


• 


0.05 



10 


• * . ■ ■ * ' 





0 

. • 

f ■ ■ " 


0 

X ■ • .. 


123456789 10 123456789 10 


n n 


Figure 2. Eigenfrequencies (left panel) and amplitudes (right panel) of the modes present in 
Figure[l] Blue circles correspond to E s = 0.1, purple squares correspond to E s = 0.01. 


Before proceeding, let us consider the relevant range for E s in hydrodynamic approxi¬ 
mation of the HVBK equations. The vorticity of the superfluid is given by the circulation 
per vortex multiplied by the number of vortices per unit area, hence 


2 fi = r x 


N 

7T R 2 


From (2.8), (2.30) and (3.19) we then obtain 


E k = 


R 2 


2 L 2 N 


In ( — 
a 0 


(3.19) 


(3.20) 


The In term is of order unity or an order of magnitude more, so that for a vessel of 
aspect ratio close to unity we have E s ~ 1/N. We require N 1 for the hydrodynamic 
approximation to be valid, hence we must have E s <C 1. 

In Figure [l] we plot the response of the container for E s = 0.1 (blue curve) and 
E s = 0.01 (purple curve) for 7 —>• 00 (perfect pinning) and K = 1 in both cases. The 
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first 500 eigenmodes of (3.17) have been plotted. We find that the container oscillates 


about the centre of mass of the system in a smooth sinusoidal manner, with a period 
much greater than the rotation period. This smooth oscillation can be understood by 
looking at the modal decomposition of the solution, plotted in Figure [2] On the left hand 
side we present the eigenvalues, in increasing order of Im(s„). On the right hand side we 
plot the corresponding amplitudes of each mode. We see from the right panel that the 
dominant contribution is from the fundamental mode, with very little contribution from 
the higher modes. As E s is reduces, we find that all the power becomes concentrated in 
the fundamental mode. Therefore, for E s <C 1, the impulsive acceleration of the container 
only excites the fundamental mode of the system. 

To obtain an analytic result for the fundamental mode, we investigate the limit £ s < 1. 
The frequency of oscillation is much smaller than the rotation frequency, hence we look 
for solutions in the limit s<Cl. We then find 


a 2k 


A 


2s 3 7 s 


(2 + k'ys) [4 (1 + K ) 7 + 2s + k'ys 2 ] 


where k = ^/2/E s >> 1. Equation (3.21) has zeroes at 

1 


s± = 


fcq . 


1 ± v/l - 4 (1 + K) /c 7 2 


(3.21) 


(3.22) 


The zero at s = —2/(A 7 ) is also a zero of A, and therefore not a pole of /. Evaluating 
R(s±) we find 


/(*) = 


1 


1 + K 


1 + 


K 


s + — s_ 


(s+e^-s.e^) 


(3.23) 


We now examine the behavior of ( |3.23[ ) by considering two limiting cases. To recover 
the oscillatory solution observ ed in Figure [l] we consider the strong pinning limit 7 
[4 (1 + K ) k} 1 / 2 . From (3.22) we have 


s± = +2i 


1 + K 


and (3.23) becomes 


m = 


1 


l + K 



We therefore obtain the period of oscillation 


tp = 


27T 


(3.24) 


(3.25) 


(3.26) 


fl|Im(s±)| Q V 1 + K 

This is the principal result of this paper. A single, long time-scale mode has emerged 
for E s < 1, which has a period proportional to E s 1 ^ 4 . This is analogous to classical 
Ekman pumping, where the Ekrnan time (proportional to A -1 / 2 ) emerges for £<C 1. 
However, whereas a viscous fluid is dissipative and results in relaxation, the superfluid 
is dissipation-less and oscillates. For E s = 0.01 we find Qtp = 8.35 in agreement with 
Figure [I] 


The weak pinning limit, 7 <C [4 (1 + K) k) 1 ^ 2 , (3.22) gives 


s+ = 


2 

&7 ’ 


2(1 + AT) 7 . 


(3.27) 
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Because s_ <C s + in this limit, (3.231 becomes 


/ 


1 

1 + K 


(1 +Ke s ~ t ) . 


(3.28) 


This exponential decay is a result of the friction force of the vortex lines as they slide on 
the boundaries. When K = 0 the timescale for the relaxation is 


1 _ 1 
fi|s_| 2127 


(3.29) 


which is the result obtained in Reisenegger (1993). Because the weak pinning limit was 
studied by Reisenegger (1993), the principal focus of this paper will be on our new result 
in the perfect pinning limit. 


3.2. Ekman pumping 

To visualize the flow for E s <C 1 we can invert the other variables in this limit. We obtain 


V* = V* + V S B 
Xs =xl + x? ■ 


(3.30) 

(3.31) 


where we have separated the solution into the interior flow and boundary layer corrections 
(denoted with subscripts I and B respectively) given by, 

o s —^ _ e p s +t ' 


v! = 


l 


v = 

’ a 


1+K 

1 

1 + K 


1 - 


s+e 


S_|_ — s_ 


s + e s+t - s_e s -‘ 


s + — s_ 


e k(z-l) _ e ~k(z+ 1) 


J _ 2z (e s_t — e s+t ) 
Xs 


xf = 


k (s + — s_) 

2 (e s t — e s+t ) 


g fc(2-l) _ e -fc(2 + l) 


k (s + — S-) 

There is no boundary layer component for the pressure, which is given by 


Ps = 


1 + K 


1 - 


s_|_e s —* — s_e s+< 
s+ — S- 


(3.32) 

(3.33) 

(3.34) 

(3.35) 

(3.36) 


Equations (3.32)—(3.36) are akin to Ekman pumping in a viscous fluid (Greenspan 


& Howard 1963 Benton & Clark 1974 van Eysden & Melatos 2013), where the expo¬ 
nentially decaying term has been replaced with the oscillatory e s±t terms. The variable 
Xs is like a stream-function for a secondary flow [see (2.22)] which draws fluid radially 


inwards in the interior, delivers it to the boundary layer where it is pumped radially 
outwards. However, when there is perfect pinning, the flow is oscillatory, with period 


given by (3.26). The Ekman circulation initially proceeds in the classical manner, clock¬ 
wise in the region 0 < 2 < 1 and counter clockwise in the region — 1 < z < 0, but is 
periodically reversed. This reversal is induced by vortex array. In the boundary layer, 
the vortex array is sheared, storing potential energy in the form of vortex line tension. 
The potential energy increases until the Ekman pumping is eventually slowed, halted and 
then reversed. As there is no dissipation in the system, this energy exchange continues 
indefinitely. Because the flow obeys the Taylor-Proudman theorem to leading order, the 
interior azimuthal flow is columnar, as in classical Ekman pumping. The vortices in the 
interior remain straight and move inwards and outwards with the secondary radial flow, 
increasing and decreasing the angular momentum of the fluid. 
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Figure 3. Dimensionless angular velocity f(t) (blue curve), azimuthal angular velocity, Vj (/,) 
(purple curve) and radial velocity (divided by r) dxl /dz (red curve) and the radial vortex line 
displacement U^ (brown curve). We take E s = 0.01, I\ = 1 and 7 —» 00 . Angular momentum is 
exchanged between the azimuthal flow and the container (blue and purple curves). The radial 
flow (red curve) is n/2 radians out of phase with the azimuthal angular velocity. The vortices 
(brown curve) move inwards and outwards with the radial flow. 


In Figure [3] we plot the dimensionless angular velocity of the container, fit) (blue 
curve), angular velocity of the fluid, V/(t) (purple curve) and radial velocity of the fluid 
(divided by r), dxl/dz (red curve) and the radial vortex line displacement (brown 
curve) for E s = 0.01, K = 1 and 7 -7 00 . The top two curves (blue and purple) show that 
the angular velocity of the container and the interior fluid are sinusoidal and 7 r radians 
out of phase. Therefore angular momentum is conserved between fluid and container, as 
required. The bottom curves show that the radial component of the secondary flow and 
the radial vortex line displacement are 7 r /2 radians out of phase. As the radial flow shears 
the vortices in the boundary layer, the restoring force from the line tension increases until 
the radial flow is halted, and then reversed. At this turning point in the radial flow the 
vortex displacement is at a maximum. The vortices are only displaced radially inwards, 
because the azimuthal flow velocity is always greater than or equal to its initial value. 

Snapshots of the flow described by Figure[3]are presented in Figure[4]at the times t = 2 
(left panel) and t = 4 (right panel) The magnitude of the fluid angular velocity V s is shown 
in color, where dark is zero and light is unity. The radial vortex line displacement is 
over-plotted as thick-black lines, scaled appropriately for illustrative purposes. Contours 
for the stream-function of the secondary flow, \ s are plotted as dashed lines for values 
of ± 0.002, 0.01, 0.025, 0.05, 0.1. In the left panel (t = 2) the secondary flow is near 
a maximum, represented by the large density of contours. The flow is clockwise in the 
upper-right quadrant and the secondary flow is moving the vortex lines inwards in the 
interior. The angular velocity of the container and azimuthal flow are matched, so there 
is no boundary layer correction. In the right panel (f = 4), the secondary flow is near 
a minimum. The radial component of the secondary flow is turning from inwards to 
outwards in the interior. The vortex lines are at their maximum extension, and will 
subsequently straighten to drive the outward radial flow in the interior. The difference 
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Figure 4. Flow plots for E s = 0.01. Left: t = 2 (near secondary flow maximum), right: t — 4 
(near secondary flow maximum). The color shading represents the magntude of the azimuthal 
angular velocity, Vj (t), where blue zero and white is unity (c.f. Figure [ 3 ]). The dashed con¬ 
tours show streamlines for the secondary flow xl(t) with values: ± 0.002, 0.01, 0.025, 0.05, 0.1. 
The solid lines show the vortex line displacement Ui, plotted with exaggerated amplitude for 
visibility. 


between the angular velocities of the container and fluid is greatest, and the boundary 
layer in the azimuthal velocity is prominent. 

To make a final connection to Ekman pumping, we examine the governing equations in 
the limit E s -C 1. In the interior, we fi nd fro m (3.32) that the azimuthal velocity evolves 
on the timescale E s and from (3.34) the stream-function for the secondary flow, 

Xs, scales as El //A . Applying these scalings, we find in exactly the same manner as for 
classical Ekman pumping [see Greenspan & Howard (1963), §5] the solutions for the 
interior flow satisfy 


0 = 2V. 1 - P„ 


0 = 


dVj , M 


dt 


dz 


(3.37) 

(3.38) 


The second, (3.38) describes the geostrophic balance in the interior. In the boundary 
layer we have 


0=2 -E, 


0=2 -E„ 


dz 2 

d 2 


V B 

v s ? 

dxf 


dz 2 J dz 


The boundary conditions are given by (2.36)-(2.38), where (2.37) is replaced by 

d d x f 


- 2 /± 


E s d 2 x ; 


B 


dt dz 7 dz 2 

at z = ±1. Finally, the initial conditions are 

TA/(0) = 0, 


P s = 0 , 


(3.39) 

(3.40) 

(3.41) 


(3.42) 
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(3.43) 

(3.44) 


Equations (3.371 (3.441 are sufficient to recover the approximate solution (3.30)-(3.36). 


They are the equivalent of the boundary layer approximation derived for classic spin-up 
by Greenspan & Howard (1963). They will be used in (jdjto derive a general solution for 
the two-fluid system. 

The mechanism described above possesses many properties similar to Ekman pump¬ 
ing. In particular, boundary layers form on boundaries orthogonal to the rotation axis, 
driving a radial outflow. This flow is replenished by fluid from the interior, generat¬ 
ing a secondary circulation. However, there are important differences. Classical Ekman 
pumping is dissipative, while the mechanism described here is dissipation-less. Also, in 
classical Ekman pumping, the Coriolis force couples the azimuthal and radial velocity 
in the boundary layer, driving the secondary flow. Here, the velocity components in the 


boundary layer are uncoupled [see (3.39) and (3.40)], but it is the boundary condition 
itself (2.19) that couples the velocity components. 


Finally, (3.37) (3.41|) can be solved to obtain a boundary condition for the interior 


flow analogous to Greenspan (1968) [see (2.6.13)] that relates the differential vorticity at 
the boundary with the secondary inflow into the boundary layer, 

hm^h ) iv ' =i ' Vx (»'-*’'>), < 3 - 45 ) 

at 2 = ±1, where v 1 denotes the velocity field in the interior. This relation applies at 
the wall in cylindrical geometry and applies only to a single fluid. 


4. Helium II 

4.1. Experiments 

The response of a superfluid-filled container following an impulsive acceleration has been 


investigated in series of experiments conducted by Tsakadze & Tsakadze (1980). In §5 of 


Tsakadze & Tsakadze (1980), plexiglass cylinders are coated with a powder to increase 


vortex adhesion, and oscillatory motion of the container is reported. The experiment 
is conducted at 1.52K with a rotational frequency of 3rads -1 . The vessel is 64mm 
in diameter, 50 mm in height and 0.2 mm thick. The response of the vessel is given in 


Figure 27c of Tsakadze & Tsakadze (1980). The magnitude of the angular velocity after 
the spin up is not given, but appears to be around 3.6rads -1 . The response comprises an 
exponential decay of 0.3rads -1 over approximately 40s, followed by steady oscillation 
with a period of 40 s, decaying over a timescale of roughly 480 s. The oscillation period also 
appears to decrease with time, probably a result of non-linear effects as the experiment 
has Rossby number e ~ (3.6 — 3)/3 = 0.2. 

Tsakadze & Tsakadze |l980) also report the following effects. The oscillation period is 


unaffected by the radius of the cylinder. This was tested by inserting coaxial cylinders of 
both 22 mm and 43 mm inside the main cylinder. However, the frequency was observed to 
double when the container height was halved, which was achieved by inserting disks into 
the cylinder. This behavior is qualitatively consistent with the columnar flow of Ekman 
pumping and the solution presented in [J3] 

To apply the our theory to the Tsakadze experiments, we need to generalize the solution 
in £[3]to the two-fluid system in [J2] For the experiment described above, the dimensionless 
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parameters at 1.52 K are |Donnelly & Barenghijil998) 



B = 1.2748, 
B' = 0.2922, 
Pn = 0.1206/9 . 


(4.1) 


The third parameter, Ma, where u± is the first sound velocity, is the Mach number and 
indicates the relative importance of compressibility effects. Clearly for helium II the 
incompressibility approximation is valid. For the experiment described above, all the 


parameters of our theory are given by (4.11 except K. In experiments involving spheres 


and no pinning, Tsakadze & Tsakadze (19801 present the steady-state spin down of empty 
and filled vessels from which one can determine K ~ 0.8 (van Eysden fe MelatospOll ). 


Unfortunately, Tsakadze & Tsakadze (1980) do not provide enough data to determine K 


for this experiment, and we take K = 1. 

4.2. Two-fluid solution 

Solving the full system of equations in ([2] poses a formidable challenge. In helium II, we 
have E <Cl and E s -c 1 but we still have E s ~ E, so that neither viscosity or vortex 
tension is negligible. Also, because B , B' ~ 1, the mutual friction coupling time between 
the two fluid components is of the order of the rotation period, and the normal fluid and 


superfluid are approximately locked together over the Ekman time (Reiscncgger 1993 
van Eysden fc Melatos||2013 2014). 


To solve the two-fluid system we use a boundary layer approximation like that used in 


§5 of Greenspan & Howard (1963) and presented in 14.3 We assume that the azimuthal 
velocity is symmetric about the z = 0 plane, and hence Xn,s are antisymmetric. In the 
interior we have 


0 = 21.;' - {Vl - vj) - Pn , 




dt 


PsB' (dxi _ dxj 
p \ dz dz 


PnB' 


0 = 2 Vl + ( V n - V s) - p s . 
dV! 


0 = 


p 

,dxl 


dt dz p 
while in the upper boundary layer we have 

d 2 p s B 
' dz 2 p 

PsB ( n „ d 2 \ dx B 

2 P 


PnB (K 1 - v/) + ^ 


dxi 

dz 


dxi 

dz 


(4.2) 

(4.3) 

(4.4) 

(4.5) 


0 = —E 


dx% _( p s B' 
dz 


V L 

v n 


2 -Es 


dz 2 ) dz 


PsB' 

2 P 


2 -Es 


dz 2 


Vf 


(4.6) 
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df_ 

' dz 2 


p s B 

P 


V 1 

' n 


2 - 


PsB'\ d\ 


P 




dz 


0 = 


0 = 


p s B 

^P 

PnB 
2 P 

-o- 

PnB 


2 -E, 


(p_ 

dz 2 


K 


p s E' 
2 P 


2-E, 


d 2 \ 5 Xs z 


dz 2 J dz 
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(4.7) 


2-E, 


<9z 2 


2 P 


1 - 


PnB' 

2 P 

2-E ( 

PnB' 

2 p 


dxf _ PnB d\n 
dz p dz 


p 


dz 2 


V B - p^v b 


2-E, 


d 2 \ d X f , PnB' dxi 


dz 2 ) dz 


P 


dz 


(4.8) 


(4.9) 


These equations reduce to (5.3)-(5.10) in Greenspan & Howard (1963) and (3.371 (3.401 
in the present paper when B = B' = 0. 


The governing equation for the motion of the container is (2.31). The boundary con¬ 


ditions for the superfluid are (2.36), (3.41) and (2.38). For the normal fluid, the required 


boundary conditions are Greenspan & Howard (1963), §5] 

Vn-f = 0 , 

dy B 
Xn = 0 
dz 

Xn=0. 


The initial conditions are (3.42) (3.44) for the superfluid and 

^( 0 ) = 0 , 


(4.10) 

(4.11) 

(4.12) 

(4.13) 


for the normal fluid. The general solution is presented in ?|A] Because of its complexity, 
we do not present the final result in algebraic form. For the numbers (4.1) and taking 
K = 1, 7 —> oo and Q = 3rads -1 , we obtain 

f(t) = 0.500000 + 0.000023e -3 - 829O8t 

[0.499977 cos(0.58827t) + 0.001126 sin(0.58827f)] , (4.14) 


+ e 


- 0 . 004182 * 


where t is in seconds. The first term is the steady-state term, which corresponds to 
the centre-of-mass of the system at 1/(1 + K). The second term corresponds to a rapid 
damping and has low amplitude. This term represents the strong mutual friction between 
the two components. The third term is an oscillation, weakly damped by viscosity. It has 
a period of 10.7s and a damping time 239.1 s. 


The result (4.14) suggests that viscosity and mutual friction have little effect on the 
final result. This is because the Ekman time and secondary flow scale as E 1 / 2 in clas¬ 
sical Ekman pumping, but scale as E |^ 4 for the oscillatory superfluid Ekman pumping 


mechanism identified in this paper. For the numbers in (4.1), we find 

E 1 / 4 = 0.088, E 1/2 = 0.006. 


(4.15) 


Therefore the superfluid mechanism is dominant. Neglecting viscosity (E = 0) and as¬ 
suming strong mutual friction coupling (E,E' ~ 1 ), an approximate solution can be 
obtained. The details are presented in (JbJ The result has the form (3.23), with s± now 
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Figure 5 . Com paris on of the general solution [(4.141, thick black curve] and the approximate 
solution [(3.231 and (B 121, thin gray curve] assuming zero viscosity (E = 0) and strong mutual 
friction couping ( B , TT ~ 1). Viscosity has a negligble effect on the period, which is the same 
for both solutions. The experimentally observed period decays from an initial maximum of 40 s. 


given by 


s± = 


1 


1 ± i /1 — 


4(1 + K ) k"/ 2 p s 


(4.16) 


The oscillatory Ekman pumping mechanism identified in [J3]is operating in the superfluid, 
with the normal fluid locked to it as a result of the strong mutual friction coupling. This 
mass loading of the superfluid results in a factor of p s /p in the spin-up time. Only the 
leading order normal fluid velocity is involved in the oscillation; there is no secondary 


flow in the normal fluid. The reader is referred to ] B.2 for the solution for all variables 
in this limit. 


A comparison between the general solution (4.14) and the approximate solution [(3.23) 
and (B 12)] assuming E = 0 and B,B' ~ 1 is shown in Figure [5] The period of the 
approximate solution is identical to the general result, suggesting it depends weakly on 
viscosity. The approximate solution is dissipation-less, whereas the general solution is 
damped on a timescale of a few hundred seconds. When 7 —► 00 , the period of oscillation 
at T= 1.52 K is 


t P = 10.7s 


3 rads 

n 


-1 


-1 


5.95 x 10- 
E„ 


—5 \ - 1/4 


0.8794 \ 1/2 
Ps/P ) 


l + K 


- 1/2 


(4.17) 


This is shorter than the maximum observed period of oscillation in the Tzakadze exper¬ 
iments, which decays from an initial maximum of approximately 40 s. The period can be 
lengthened by decreasing the frictional pinning force (decreasing 7 ), however, this quickly 
shortens the damping time to much less than that observed in the experiments. Because 


of the precision to which the parameters in (4.1) have been measured, the discrepancy 


must arise either because the linear hydrodynamic theory is invalid, or one or more of the 
experimental parameters have been erroneously reported. The former is certainly likely; 
as discussed in §4.1[ the Rossby number is 0.2, so non-linear effects may be present. The 
observed decay of the oscillation period cannot be captured by the linear theory. The im- 
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pulsive spin-up of the container cannot be done instantaneously as required for a perfect 
comparison with theory, and this may explain this initial rapid decay of the container 
motion observed in the experiments. Errors in reporting the experimental parameters 
is also possible, because Tsakadze & Tsakadze (1980) report some inconsistently (e.g., 
vessel radii), and others not at all (e.g., the vessel moment of inertia, which we have 
assumed equal to that of the fluid). A thorough comparison between theory and exper¬ 
iment will require new experiments. Such experiments could be combined with the well 
known properties of helium II to provide a useful test of both the HVBK theory and 
the calculations in this paper, as in van Eysden & Melatos (2011). Alternatively, the 
container response predicted here may be a useful check of known superfluid parameters 
in helium II, or used to constrain parameters in condensates where they remain unknown 
( van Eysden fe Melatos| 2012). 


4.3. Tkachenko oscillations 


In our calculations, we have solved the complete equations of Chandler & Baym (1986), 


in order to apply the most sophisticated hydrodynamic theory available to the helium II 


experiments. Although this theory includes the extra term (2.9) over the HVBK theory, 


we find that the contribution from this force vanishes. This result was also obtained 


by Reisenegger (1993). The reason for this lies in the ansatz (2.21), where the flow is 


assumed to be axisymmetric and the radial dependence of azimuthal velocity matches 
the boundary conditions. The equation of motion for the vortex lines requires that £ has 


the same radial dependence as the superfluid velocity, and upon substituting (2.21) into 


(2.9), we find that er = 0. 


Physically, (2.9) describes the restoring force experienced by the vortex array as a 


result of displacements of vortices from their equilibrium configuration. Hence, the vor¬ 
tex array experiences a restoring force in response to shear deformations, giving rise to 
Tkachenko oscillations. However, in the spin-up problem considered here, the top and 
bottom boundaries are rigid bodies and the impulsive acceleration of the boundary does 
not displace any vortices from their equilibrium configuration. Therefore we do not neces¬ 
sarily expect Tkachenko modes to be excited in the experiments ofjTsakadze & Tsakadze| 

(fiM)|). 


Tkachenko modes may be excited at the cylinder wall, which has been neglected in our 
analysis. However the required boundary conditions required at the cylinder wall are not 


discussed by Baym & Chandler (1983). It is known that it is energetically favorable for 


vortex-free region to form adjacent to the cylinder wall, making it difficult for changes 
in the angular velocity of the side-wall to be communicated to vortex array. Even if they 
were excited, Tkachenko oscillations on timescales of order much less than the 

Ekman pumping mechanism described above, which has timescale E s . Therefore 

the motion of the vortex lines is determined by the vortex tension producing Ekman 
pumping, and the restoring force from the lattice deformations is too slow to have any 
effect. 

Therefore, Tkachenko oscillations are unlikely to be excited in the experiments of 


Tsakadze & Tsakadze (1980). Tkachenko oscillations have been imaged in rapidly-rotating 


dilute-gas Bose-Einstein condensates, excited by blasting atoms or creating a dip in the 
trapping potential at the centre of the condensate (Coddi ngton et aL|2003 ). Here, the ex¬ 
perimental setup and excitation mechanism are appropriate to excite Tkachenko modes, 


whereas in the experiments of Tsakadze & Tsakadze (1980) they are not. Therefore, the 


Ekman pumping mechanism presented here provides a much more natural explanation 


for the oscillations observed in the Tsakadze & Tsakadze (1980) experiments. 
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5. Neutron stars 


C. A. van Eysden 


Neutron stars are compact stellar corpses formed by the core collapse of a massive star 
after a supernova. The resulting compact object typically has a radius of 12-14 km and 
a mass of 1.4-2 times that of the sun. The outer kilometer of the star is a rigid, highly 
conducting crustal lattice interspersed with a neutron superfluid, while the core comprises 
approximately 95% superfluid neutrons, 5% electrons and superconducting protons, and 
a small fraction of muons. In the deep core, the composition remains uncertain. 

Neutron stars are also endowed with strong magnetic fields, ranging from 10 8 -10 15 
Gauss. In isolated radio pulsars, radio waves are beamed along the magnetic axis and 
create a lighthouse effect that is observed as pulsations on Earth. Occasionally, these 
objects are observed to ‘glitch’, where the r otati onal frequency of the star suddenly 
increases by up to one part in 10 5 (Espinoza et al. 2011). The distribution of glitch 


sizes and waiting times is consistent with the hypothesis that they arise from a self- 


organised, collective, avalanche process (Anderson & Itoh 1975 Melatos et al. 2008 
Warszawski & Melatos 20131. During a glitch, angular momentum stored in the decoupled 


superfluid components of the neutron star is transferred erratically to the solid crust 


(~ 1% of the total moment of inertia) (Link et al. 19921, although the exact source of 
this angular momentum is uncertain ( Andersson et «!.||2012[ Chamel||2013 ). The trigger 
is often ascribed to collective unpinning of quantized superfluid vortices, induced when 
crust-superfluid differential rotation (and hence the Magnus force) exceeds a threshold 
( Baym et aL||1992 Warszawski et aL||2012 ). 

The glitch event is followed by a recovery phase, during which the components inside 


the star adjust to a new equilibrium configuration (Baym et al. 19696). This recovery 


is typically dominated by a step increase in the rotational frequency of the star, a step 


change in the frequency derivative, and a quasi-exponential relaxation (van Eysden & 
Melatos|'2010). However, quasi-exponential oscillations have also been reported following 


some glitches. A post-glitch oscillation with a period of approximately four months was 
reported in the Crab pulsar soon after its discovery (Rudermanjl970ai6). The 1975 and 
1986 Crab glitches were observed to “overshoot” their final rotation frequency during 
their recovery (Wong et al. 2001 van Eysden &; Melatos|2010 ), which could be potentially 
be explained by a damped oscillation with a period of approximately 200-300 days. 
Oscillations also appear to be present in timing residuals following the 1988 Christmas 
Vela glitch [see figure 2 of McCulloch et al. 1990)], which, if real, appear to have a 
period of ~ 20 days. These oscillations have received little attention observationally or 
theoretically, but may shed light on the interior of neutron stars. Here we consider the 
possibility that the superfluid oscillations described in the previous sections may be 
operating in neutron stars, and ask what the observable consequences are for glitch 
recovery. 

During glitch recovery, the proton-electron plasma in the core is coupled to the highly 


conducting crust via the magnetic field on a time-scale of seconds (Easson 1979 


van 


Eysden 20141, and hence, these components are assumed to be rigidly locked together. 


The superfluid neutrons are usually assumed to respond changes in the angular velocity 
of the crust via the interaction of the neutron vortices with the proton-electron plasma 
in the core ( Alpar et a7||1984 Mendell||1991a|6 Glampedakis et alj|201 1 1 . The electrons 
scatter from the neutron vortices that are magnetized by the entrainment of proton 
currents, giving the mutual friction force 


F = -—/3 0 lo s x [w s x [v L - n„)] , 
Pn 


(5.1) 


where the subscript n refers the the proton-electron plasma, which plays the role of the 
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normal fluid component and s refers to the neutron superfluid. The coefficient /3 q is given 


by Mendell 1991 &) 


2 (y-i ) 2 


, 7/6 


A) = l-lx 10- 2 i /2 , , , (5-2) 

y L ' z (1 — x) 

where x is the proton fraction (p n /p in the notation used here) and y is the normalized 
effective mass of the proton from Fermi liquid theory. The coefficient f3' 0 that appears 
in (5.1l in HVBK theory is typically taken as zero in neutron stars. The motion of 


magnetized neutron vortices is also inhibited by their pinning to flux tubes. When a 
vortex is pinned it moves with the normal fluid so that = v n . However, an unpinned 


vortex moves with velocity given by (2.11). Solving the vortex line equation of motion 


(2.11) gives 


VL = V n - 


1 


{w s x [w s x (v s - v n ) + <y + 1\ 


1 + Po 

+ Po [w s x (v s -v n ) + <r + t]} . 


(5.3) 


At zero temperature, all vortices are pinned to flux tubes. At finite temperature a vortex 
can become thermally excited with activation energy A upon unpinning. The activation 
energy for unpinning depends on the magnetic energy between a vortex and a flux tube 
E p , the dimensionless vortex tension 7~, the angular velocity lag Aw and the critical 
angular velocity lag for unpinning Aw c as 

A = 5.LE P T 1/2 (l-p- 
V Aw, 


5/4 


(5.4) 


Therefore the mutual friction force has a non-linear dependence on the velocity through 
Aw. However, to simplify our analysis we take A as constant, noting that a rigorous 
calculation should include this non-linear dependence. 

The partition function for this two-state system is 

o-P A 


Z = 1 


(5.5) 


where /? = (fc^T) 1 , ks is Boltzmann’s constant and T is the temperature. When the 
thermal energy is much greater than the activation energy for unpinning a vortex /3 _1 


A , the all vortices are thermally activated and move with average velocity given by (5.3). 


However, when the thermal energy is much less than the activation energy for unpinning 
/? _1 <C A, the vortex movement is exponentially suppressed pinning the vortices to the 
normal fluid. For slow vortex slippage, e~^ A 1 and Z sa 1 and the ensemble average 
vortex line velocity is given by 


(vl)i 3 ~v n = {v L - v n )p 

1 

{w., 
J o 


l + £o 2 


« z : e 0A (v L - v n ) 
x [w s x (v s - v n ) + er + t] 


+ p 0 [w s x (t> s -v n ) + cr + t}} e 


-PA 


(5.6) 


Substituting the vortex line velocity from (|5.6|) into the vortex line conservation equation 


(2.10) and integrating, we obtain an HVBK-likc equation for the superfluid (2.2) with 


modified mutual friction coefficients 


B = 


B’ = 


2P Po —PA 
Pnl+Pl 
2 p 
Pn 


1 - 


i + 


o-PA 


(5.7) 


(5.8) 
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For a detailed derivation of the theory of vortex slippage, the reader is referred to |Link| 
(2014), where (5.6) and (5.8) are rigorously derived [see (48) and (51) in that reference]. 

Before proceeding, we assess the validity of the incompressibility assumption for the 
present application. In neutron stars the speed of sound is approximately one fifth of 


the speed of light (see e.g., Reisenegger & Goldreich (1992)), giving a Mach number 
1.7 x 10 -3 . Therefore the incompressibility assumption in neutron stars is valid. 

We can now apply the theory in (j 2 ]to neutron stars. Because the normal fluid is rigidly 
locked to the crust via the magnetic field, we have v n = tl c r. Therefore, it is a single- 
fluid problem involving only the superfluid and we can repeat the analysis in [JdJwith the 
addition of mutual friction terms. In this case, the crust plays the role of the container. 
We assume that the vortices are perfectly pinned to the crust with 7 —>• 00 at 2 = ±1. 
In the core of the star (the region — 1 < z < 1), we expect superfluid oscillations as 
illustrated in m where vortex lines are straight in the core, except a small boundary 
layer region adjacent to the crust. However, the motion of the neutron superfluid in the 
core is now impeded by the interaction with the flux-tubes in the proton electron plasma, 
which is described by the mutual friction force. 

For the superfluid oscillation mechanism to operate, vortices pinned in the crust must 
communicate with vortices in the core via vortex line tension. In the outer crust, neutrons 
pair in a 1 So state, while in the core they are expected to pair in a 3 P 2 state. The 1 Sq 
phase is expected to be lost near the base of the crust, however current calculations of the 
pairing gaps for the 3 P -2 phase are not certain enough to know at what depth superfluidity 
will be restored. If the two phases are connected, then the relevant boundary condition 
at the interface are also unknown. It may be that the vortices in the crust and core are 
not connected, in which case the superfluid oscillations described here will not occur. 

Repeating the analysis in ^3] with the mutual friction terms, the result (3.21) generalizes 
to 


ke 


2k 


2 ps (ps + p„.B) 

1 

and the inverse Laplace transform is 
1 


{(2 p-p n B') [2Kp s + (2p-p n B') 
(ps + p n B) (pks + p n B)} , 


m = 


where 


(2 p - p n B') + 2Kp s 


(2 P — PnB') + 


2Kp s 

S+ — S- 


(s+e s -‘ - s_e s+t ) 


s± = 


PnB k ( PnB ) 2 - 4 (2p - p n B’) [2Kp s + (2 p - p n B')\ 


2 P 


± 


4 p 2 k 


(5.9) 


, (5.10) 


(5.11) 


In a typical neutron star, we have p n /P = x ~ 0.05, f l = 100rads -1 , 0.3 ^ y = m*/m p ^ 
0.7 and K ps 50 (Mendell 1991 a | 6 |van Eysden fc Melatos 2010). The pinning fraction 
in the core is uncertain. Hence we have (3q <C 1 [see (5.2)], and also Kp s "^> p because 
the neutron superfluid is the dominant contribution to the moment of inertia of the star. 
The superfluid Ekman number is approximately 


E s = 3.06 x 10 


-17 


( - - -) 

\ 100 rads -1 / \ 10 6 cm 


(5.12) 


Under these assumptions (5.10) can be written 

27rf' 


fit) = e 




cos 


tr) 


2nt d 


sm 


2irt 


(5.13) 






































Figure 6. Damping time (left) and period (right) for oscillations in a neutron star. Lines are 
plotted for y = m*/m p = 0.3 (blue) and 0.7 (purple). In the right panel the period is independent 
of rrip and the lines overlap. 


where 


1 

^ d fi|Re(a±)| 

o eP A 


(5.14) 

27 T 

Tre ^ A / 2 

1 pk 

(5.15) 

p f2|Im(s±)| 

n \ 

I<Ps' 


In Figure[6]we plot the damping time and period for superfluid oscillations in a neutron 
star as a function of the pinning parameter (3 A. We plot the full result for the damping 
time and period using (5.11). Results are plotted for y = m* /m p = 0.3 (blue) and 0.7 
(purple). We find that the approximate results given in (5.14) and (5.151 give excellent 
agreement. The damping time becomes longer as pinning is increased, because the amount 
of vortices moving and dissipating energy through scattering is exponentially suppressed. 
The period also increases as pinning increases, because pinning is inhibiting the movement 
of vortices and hence the transport of angular momentum. In the limit of infinitely strong 
pinning, the vortices stop moving and the period becomes infinitely long. 

In principle, oscillations of almost any period are possible, depending on the pinning 
strength, A. The oscillation amplitude is unaffected by pinning and always equal to 
the glitch amplitude. The damping time-scale is always longer than the period, so the 
oscillations are only weakly damped and the second term in (5.13) is small. Post-glitch 
oscillations are observed in both the Crab and Vela pulsars that are candidates for the 
oscillatory Ekman pumping mechanism in this paper. In the Crab, oscillations with a 


period of roughly four months were reported following the earliest glitches (Ruderman 
1970a b). The “overshoot” observed in the recovery the 1975 and 1986 Crab glitches could 


also correspond to a damped oscillation with period of 200-300 days ( Wong et q/.||2001 
van Eysden & Melatos 2010). The former would require /3 A « 23.8, and for the latter 


/3 A « 24.8. In Vela, a damped periodic oscillation is clearly visible in the timing residuals 
following the 1988 glitch [see Figure 2 of McCulloch et al. (1990)]. These oscillations 
appear to have a period of ~ 20 days, requiring (3A ~ 20.2. Therefore, superfluid Ekman 
pumping can explain these oscillations if pinning between vortices and flux tubes is 
extremely strong, limiting vortices to a very slow creep. 

These numbers can be compared with theoretical expectations using known parame¬ 
ters in the outer core of the pulsar. Assuming a typical pulsar temperature and using 
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parameter estimates from Link (20141, from (5.4) we find 


PA = 2.05 x 10 


5 , 100MeVV / V0- 12 \ 1/ V 1 ° 7KS 1 


E n 


T 


T 


(5.16) 


where we have taken A u = 0. Clearly (5.16) leads to extremely long oscillation periods 


and damping times. However, this is a gross upper estimate because the lag Aw is not 
likely to be zero. If the lag is close to the critical lag for unpinning A« c , the estimate 


(5.16) can be greatly reduced, see (5.8). However, to correctly determine the lag the 


non-linear calculation including velocity dependence in the mutual friction coefficients is 
required, which is beyond the scope of this paper. 

Finally, we remark on limitations of the present model. Neutron stars are strongly 


stratified, typically with buoyancy frequencies N exceeding the rotation rate (Reiseneg- 


ger & Goldreich 1992). Buoyancy inhibits Ekman pumping from penetrating the core, 
confining the circulation to a buoyancy layer with thickness flL/N (Walin 1969; Abney 
fe Epstein|1996[|van Eysden fe Melatos|2008 ). The spin-up is time is reduced commensu- 
rately reduced to E~ 1/ ' 1 N~ 1 because the effective size of the Ekman cell has been reduced 
by a fraction Cl/N. The final state is not co-rotation, but has persistent shear between the 
fluid in the buoyancy layer and that in the core which in untouched by Ekman pumping 
( Melatos|2012 ). In a typical neutron star, we find Cl/N 0.20 (Reisenegger & Goldreich 
1992), so that the spin-up time is reduced by a factor of 5. 


A similar result is expected for the superfluid Ekman pumping mechanism discussed 
in this paper. As in the classical stratified Ekman pumping, we expect that buoyancy 
forces will confine the circulation to a depth CIL/N. The Taylor-Proudman theorem no 
longer applies in the core, and the vortex lines will be sheared. Above CIL/N the fluid 
undergoes in superfluid oscillations, but at greater depths the fluid does not participate. 
To estimate the period of oscillations, we assume that the length scale for the Ekman 


cell is reduced from L to CIL/N , and (3.26) becomes 


t P — 


7T Cl 

n \l n \ i + k 


(5.17) 


In this case we have yT T/N ~ 0.45, shortening the period by a factor of approximately of 
2.2. Therefore the superfluid oscillations are more weakly affected by stratification than 
classical Ekman pumping. 


6. Conclusions 

The response of a uniformly rotating superfluid threaded with a large density of vor¬ 
tex lines to an impulsive response of its container has been investigated. When vortices 
are strongly pinned to the container and there is no dissipation, the system oscillates 
persistently with a period of order where E s is the dimensionless vortex line 

tension parameter. The low-frequency mode is generated by a secondary flow analogous 
to classical Ekman pumping, periodically reversed by the vortex tension. This secondary 
flow transports vortices radially inwards and outwards in the interior, changing the an¬ 


gular momentum of the fluid. The full set of equations of Chandler & Baym (1986) have 


been solved and compared with experiments performed by Tsakadze & Tsakadze (1980) 


in superfluid helium. We find qualitative agreement between theory and experiment, 
however quantitative comparison is difficult because some experimental parameters are 
unknown, and the experiment may probe the non-linear regime, which is beyond the 
scope of this theory. The columnar nature of superfluid Ekman pumping also reproduces 
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the experimentally reported dependence of the frequency on vessel radius and height, 
however the theory would benefit from comparison with new, carefully conducted exper¬ 
iments. We argue that Tkachenko modes are unlikely to be excited by the boundary and 
initial conditions imposed in the spin up experiments of Tsakadze & Tsakadze (1980). 
In neutron stars, the movement of neutron vortices is inhibited by pinning to type II 
superconducting protons in the core, slowing the superfluid Ekman pumping process. 
For very strong pinning, and hence slow vortex creep, the theory can explain post-glitch 
oscillations with periods of days to years, with much longer damping times. 


The author thanks Bennett Link for useful feedback on the manuscript. 


Appendix A. General solution in the boundary layer approximation 

Here we present the solution to the governing equations (4.2)—(4.9). The solution is too 
large to present algebraically, so we outline the method and present key results. Because 
if the size of the equations, it is easier to perform the calculations in programs such as 
Mathematica or Maple and handle them numerically. 

We assume the container is rigid so that the upper and lower boundaries are move 
identically. Therefore, the V U:S are symmetric in z and \n,s are anti-symmetric. Applying 
these boundary conditions, we solve in the upper-half plane only. 


The normal fluid has the solution to the interior flow equations (4.2) -(4.5) 

3 


Vn(t) = J 2 V r 


1 e Sit 


E 


n4 > 


(Al) 


and similarly for V s . Three time-scales exist in the problem and are obtained when 
applying the boundary conditions below. The stream-function for the secondary flow is 


xh{z,t) = x 


and similarly for \s- The pressure is 




' Xn4 


(A 2) 


P n {t) = Y J Pme Sit + Pr 


n4 > 


(A3) 


and similarly for P s . Substituting (A1)- (A 3) into (4.2)-(4.51 allows to solve for the 
coefficients of Xn,s and P n . s in terms of those of V r [ s . 


In the boundary layer, the solutions to (4.6) (4.9) have the form 
3 


Vn ( z >t) = 

2=1 

V s B (z,t) = J2 


V£+e" Ml " z) + ^ l_e- k - {1 ~ z) 


V B ' 

v S2 + 


-fc+(l-*) + V B e -k.(l-z) + yB e -k(l-z 


(A 4) 
(A 5) 


in the upper boundary layer, where 


k± 


±2 ip [B±i( 2- B')\ 


E [p n B ±i(2p- p n B')\ 


k = 



(A 6) 
(A 7) 
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For Xn s we have 


XnOM) = X/ 


4=1 

3 


V -B p-fc+C 1 -^) _L -fc_(l-2) 

Ani+ C ' ' A ni—*-' 


xf(z,t) =J2 eSit xf l+ e ~ k+{1 - z) +X*-e- k - (1 ~ z) +xf,e- fe(1 - 2) 


(A 8) 
(A 9) 


Note that for the classical Ekman boundary layers (with exponents fc±), the superfluid 
and normal fluid components are coupled, whereas, for the superfluid Ekman pumping 
mechanism, the normal fluid is not involved. Substituting (A4) (A9) into (4.6 1 (4.91 
and comparing like terms gives xf± an d xfi_ 

Vgf and xfi are unrelated to other boundary layer coefficients. 

The motion of the container has the solution 


in terms of Vj^ ± 


The coefficients 


f(t) = f* 


A, 


(A 10) 


where the fi can be determined in terms of V^ ± , Vj? and by substituting (A4)— 


(A9) into (2.31). Note that only the boundary layer corrections to the velocity contribute 


to the torque on the container. 

The coefficients V^ i± , V s f and xfi an( l Vfi an d V fd can be sol ved fo r in terms of 
V ] nl through application of the boundary conditions (2.36), (3.41) and (2.38) for the 


superfluid and (4.10) (4.12) for the normal fluid and comparing like terms. For each i, 


the six boundary conditions determine five unknowns in terms of V.i , the final equation 
is an eigenvalue equation for s*. The equation is cubic is s t giving solutions for the three 
Si. The boundary conditions also give /,i = V S 4 = V n 4 . 

The remaining four unknown coefficients V)! and V n 4 are determined by the initial 


conditions (3.42) (3.44) and (4.131. The result is too cumbersome to present algebraically, 


but is presented for the experimentally measured parameters in helium II in equation 


(4.14). 


Appendix B. Solutions in the boundary layer approximation for 
strong coupling 

In helium II, we have B,B' ~ 1, resulting in the strong coupling between the two 
fluid components. Any differential rotation is removed over the rotational time-scale 
and the two fluid-components are “locked together” over the much longer Ekman time 
(Reiseneggcr 1993 van Eysden & Mclatos 20131. This assumption can be used to find a 
more analytically tractable solution to that presented in jjA| 


When B = B' = 0, (4.3) and (4.5) describe the geostrophic balance in the interior, 


balancing the time rate of change of the a zimutha l flo w with secon dary flow. In (4.3) 
all terms are of order E 1 / 2 (Greenspan & H oward| jl963) and in (4.5) all terms are order 
(see g2|). However, when B,B' ~ 1 the third terms in (4.3) and (4.5) are order 


unity, while the remainder are of the order of the secondary flow. Therefore, we perturb 
the leading order velocity so that to leading order we have Vf = Vj and the velocity 
difference SV 1 is of the order of the secondary flow, i.e., 


V = Vi - SV 1 . 


(Bl) 
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Equations (4.21 (4.51 become 


0 = 2U„ J - P n , 
n _ SVn . pdxj 
dt + “ dz 


P P \dz dz) ’ 


0 = 2E/ - P s , 

n ^V/ , o d xl PnB XJrI , p n D , (dXn \ 

0 = ^T + 2 a7 _ ^ sy + 7 


The solution for has the same form (All, and similarly for 6V 1 we have 


SV^^SV/e^ + SVi, 
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(B2) 
(B 3) 

(B 4) 
(B5) 


(B6) 


Using the same form for the stream-functions (A2) and p ressu res (A31, their coefficients 
can be determined in terms of V ™ and SV/ using (B 2) (B5). 

The remaining procedure follows as in ( [A] however we find that there are now only two 
time-scales si and S 2 ', the third corresponded to mutual friction coupling which is now 
effectively instantaneous. Accordingly, we can no longer impose separate initial conditions 
on Vj and V^, which must be the same. The resulting solutions are still algebraically 
cumbersome, so we present results in the limit of negligible viscosity and negligible vortex 
tension. 


B.l. Solution for negligible viscosity 

A useful result discussed in §4.2| is the case when viscous effects are negligible. We argue 
in that section that this is relevant in Helium II. 


In the boundary layer we find that V^ si ± 


= 0 and Xn,si± = 0- We also do not need to 


satisfy the normal fluid boundary conditions for no slip (4.10) and (4.12). The result is 


where 
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The implications of this result are discussed in §4.2| 
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B.2. Solution for negligible vortex pinning or vortex tension 
The focus of this paper has been on pinned vortices, however, as a useful check we can 


recover the previous results of Reisenegger (1993) and van Eysden & Melatos (2013 


2014| where pinning is negligible. 

When 7 —> 0, we obtain the result 

/(*) = TTk ( 1 + KeS3t ) ’ 
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V 1 = _-_(l - e S3 ‘) 

n 1+K V >’ 

pe S3t (kf + fcf) [(fci + k 2 ) + i/3 (fci - k 2 )\ 


5V 1 = 


4 /3 2 klk%p s 


j _ ize S3t (fci - k 2 ) 
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X s Xn 


where 


2kik,2 
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A/3 2 klklp s 
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S 3 — — — (fci + k 2 ) (1 + I\) p n 


(B 14) 
(B 15) 
(B 16) 
(B 17) 
(B 18) 

(B 19) 


Interestingly the same result is obtained by taking E s = 0. This is just a statement that 
when there is no vortex pinning, the vortex tension plays no active role in the spin-up, as 


discussed in van Eysden & Melatos (2013). There is a slight difference in the boundary 
layer solution; for E s = 0, Vjf = Xs 3 = 6 : while for 7 = 0 they are non-zero. 


From (B 19) we obtain the spin-up time 
1 P 


t s = 


fi|s 3 | PnEW/l+K) 


B 2 + (2-B'f 
(p n B) 2 + (2p- p n B'f 


1 1/2 


- 1/2 


2 p s pB 


(,PnBf + (2p - p n B'f 


(B 20) 


This result was obtained by Reisenegger| (19931 and explored in detail by van Eysden 
& Melatos (2013). The term in the curly braces arises from the coupling between the 


two fluid components is a boundary layer of order unity. The pre-factor is the familiar 
Ekman time modified by the normal fluid density fraction. The two fluid components are 
strongly coupled (i.e., they achieve co-rotation within a rotation period as B ~ 1), but it 
is only the normal fluid providing the Ekman pumping, which is reflected in the spin-up 
time. 
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